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SECTION  1 
INTRODUCTION 


i 


One  important  part  of  the  total  external  SGEMP  problem  is  the 
calculation  of  structural  response.  Since  the  structural  response  can 
couple  into  critical  subsystems  it  is  imperative  to  know  what  frequencies 
to  expect  from  structural  resonance  and  how  long  an  electrical,  structural 
disturbance  will  persist.  This  report  addresses,  primarily,  the  question 
of  radiation  damping  for  a particular  type  of  structure.  A general  point  of 
view  is  taken  throughout  the  report,  however,  so  that  the  methods  used  for 
finding  frequencies  and  damping  constants  can  be  applied  to  the  broader 
class  of  systems  to  which  the  particular  system  analyzed  belongs.  The 
particular  system  analyzed  is  similar  to  a system  like  the  FLTSATCOM  satel- 
lite or  the  RESMOD  model  (see  Figure  1)  used  in  the  MRC  Phase  IV  A and  B 
exploding  wire  experiments. 

Damping  depends  upon  frequency.  Calculating  the  frequencies  for 
RESMOD-type  systems  has  been  addressed  by  Marin1  and  Higgins2  utilizing 
lumped  parameter  circuit  models.  In  Section  2 we  set  up  the  equations  neces- 
sary to  find  the  resonant  frequency  of  N capacitive  bodies  and  then  particular- 
ize the  equations  to  the  three  capacitive  body  RESMOD  system.  The  actual 
RESMOD  structure  is  used  as  a numerical  example.  Since  the  damping  rate  is 
related  to  the  time  rate  of  change  of  energy  of  the  system,  we  also  set  up 
the  equations  describing  the  energy  of  an  N body  system.  Again  we  particularize 
the  equation  to  the  appropriate  three  body  system  using  simplifying  approxi- 
mations. The  N body  system  is  chosen  for  discussion  because  it  reemphasizes 3 
the  more  general  applicability  and  usefulness  of  the  method.  It  is  also  a 
vehicle  for  pointing  to  the  assumptions  of  this  type  of  analysis. 
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Approximations  and  calculations  for  the  damping  constant  of  the 
lowest  two  resonant  modes  of  the  RESMOD-type  system  are  made  in  Section  3. 

For  systems  which  do  not  damp  too  appreciably  in  one  cycle  the  methods  used 
here  are  quite  general.  Essentially  the  method  consists  of  making  an 
educated  guess  as  to  what  the  current  distribution  on  the  various  parts  of 
the  system  are  and  then  integrating  the  related  Poynting  vector  over  a sphere 
at  infinity.  Graphs  resulting  from  the  calculations  are  also  presented  in 
Section  3.  These  graphs  can  be  used  to  calculate  the  damping  of  symmetric 
and  asymmetric  modes  of  general  RESMOD  or  FLTSATCOM  like  systems. 
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SECTION  2 

EQUATIONS  DESCRIBING  LUMPED  PARAMETER  RESONANT  SYSTEMS 


In  this  section  we  will  first  set  up  the  equations  describing 
the  resonant  frequencies  of  N capacitive  bodies  connected  by  N-l  inductive 
bodies.  The  capacitive  bodies  are  associated  with  all  the  electric  field 
energy  and  the  inductive  bodies  are  associated  with  all  the  magnetic  field 
energy.  This  approximation  to  an  actual  structure  is  valid  if  the  distri- 
buted capacitance  of  the  inductive  bodies  is  small  compared  to  the 
capacitance  of  the  capacitive  bodies  and  if  the  distributed  inductance  of 
the  capacitive  bodies  is  small  compared  to  the  inductance  of  the  inductive 
bodies.  The  wavelength  of  the  resonant  frequencies  found  by  this  method 
must  be  larger  than  twice  the  geometric  length.  In  Section  2.2  we 
particularize  these  equations  to  a system  of  three  capacitive  bodies  con- 
nected by  two  inductive  bodies.  This  latter  system  is  close  to  the  system 
we  actually  wish  to  analyze.  Approximations  are  then  made  which  correspond 
to  other  treatments  of  the  system.  In  Section  2.3  we  approximate  the  fre- 
quencies of  the  actual  RESMOD  and  in  Section  2.4  we  express  the  energy 
associated  with  the  N body  system,  particularizing  to  the  RESMOD  approximati 

2.1  N BODY  SYSTEM 

For  an  N body  system  (see  Figure  2)  Greei.'s  theorem  indicates  that 
N 

vi  = E Ai  Q , i = 1,N  , (2.1) 

j = l J J 

where  V.  is  the  potential  of  the  i**1  body,  is  the  charge  on  the  i**1  body 


and  A. ■ are  constants  (elastances)  dependent  upon  the  geometry.  Equation 
2.1,  a static  field  approximation,  is  valid  only  in  the  limit  of  long 
wavelengths.  The  time  rate  of  change  of  voltage  difference  between  any  two 

*.  U 4-  U • 

is  obtained  from  (2.1)  as 


bodies  (the  itkl  and  the  k^1) , AV.^ 


N 

A E 

j = l 


(A.  .-A.  ,)Q. 


L., 

ik 


d2lik 

,2  2 

d t 


i t k , 


(2.2) 


where 


AV..  = L.. 

ik  ik 


d l.. 
ik 
~> 

dt~ 


(2.3) 


and  L..  is  the  inductance  of  the  link  between  bodies  i and  k and  I.,  is  the 

lk  th  th  lk 

current  running  from  the  i to  the  k body.  We  note  that  I..  = - I,  - and 

i k ki 


L., 

ik 


A dot  above  a quantity  denotes  a time  derivative.  In  Equation  2.2 
we  are  assuming  that  there  is  no  mutual  inductance  nor  any  resistance  between 
any  of  the  bodies.  The  effect  of  radiation  resistance  on  the  currents  is  es- 
sentially treated  as  a second  order  effect.  In  order  to  find  the  resonant 
frequencies  from  Equation  2.2  we  must  relate  charge  to  current.  We  assume  the 
continuity  equation  to  hold  for  all  N bodies: 

% - £ Itj  . U.4) 

J = 1 


j^i 

where  positive  currents  are  defined  as  those  flowing  out  of  the  capacitive 
bodies.  One  constraint  results  from  the  conservation  of  total  charge  of 
the  system,  namely 


N N 


(2.5) 


The  modes  of  the  system  are  defined  by  solving  the  system  of  equa- 


tions obtained  by  assuming  an  e 
N N 


lWt 


time  dependence: 


Y'  T (A.  . - A,  • + <jj“L  .5  . <5 .. ) I . 

s^2  1J  kj  Sj  si  jk  Sj 

All  ik  pairs  i / k 


(2.6) 
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for  all  capacitive  body  pairs  ik.  Equation  2.6  for  a many  body  interaction 
was  derived  using  very  simple  concepts.  Those  concepts  and  approximations 
limiting  the  equation  were  made  explicit  in  the  derivation.  It  should  be 
noted  that  the  only  interaction  between  all  the  bodies  is  capacitive  or 
electrostatic . 

2.2  3-BODY  SYSTEM 


The  RESMOD  type  system  consists  of  three  capacitive  bodies  con- 
nected by  two  inductors  as  depicted  in  Figure  1.  Letting 


L12  ~ L1  ’ 


L23  ~ L2  ’ 


*12  = *1  ’ 


(2.7) 


and 


*23  5 *2  ’ 


ai  ~ A11  + A22  ' 2A12  ’ 


a0  A23  ‘ A13  + A12  ' A22 


(2.8) 


a2  ~ A22  + A33  " 2A23  ’ 


we  obtain  from  Equations  2.6  the  equations 


IjCa^LjU)  ) + I2(aQ)  = 0 , 


riao 


I2(-a2~L2U)  ^ = 0 > 


(2.9) 

(2.10) 


or  the  equation  defining  the  frequencies  is 


(a^L^2)  (a2-L2to2)  - a2  ■=  0 


(2.11) 
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If  the  system  is  symmetric  about  the  plane  through  body  2 then  L = = L_ 

and  ■=  a^  and  the  frequencies  are  given  by 

, ,1/2 

■ 


and 


1^  = - I.,  (symmetric  mode) 

,1/2 


(2.12) 


(2.13) 


I j *=  + I (asymmetric  mode)  . 


From  Equations  2.8,  noting  the  symmetry  we  have 
ai  ' a0  = 2A22  + All  ' 4A12  + A13 


ai  + a0  = A11  ‘ A13  • 


(2.14) 

(2.15) 


If  we  denote  the  capacitance  between  bodies  1 and  3 by  C^-.  and  between 
bodies  1 and  2 by  Cp  we  have 


t2ci3>~‘  • #n  • An  • 


<C12> 


-1 


An  + A22  ' 2A12  ’ 


(2.16) 

(2.17) 


(The  capacitance  between  two  bodies  1 and  2,  say,  is  defined  here  as  the  difference 
in  potential  divided  into  the  charge  on  body  1 when  the  charges  on  all  bodies 
except  1 and  2 have  been  set  to  zero  and  in  addition  Q,,  = - Qj.)  Then 


a . a _ _2_  _1 4C13~C1 2 

1 ° C12  ' 2C13  ' 2C12C13 


(2.18) 


ci , -f  a c — 

1 0 2C 


13 


(2.19) 


Substituting  (2.18)  and  (2.19)  into  Equations  2.11  and  2.12  we  find  that 


12 


(2.20) 

(2.21) 


In  the  limit  that  the  presence  of  one  body  does  not  influence  the  potential 
on  the  other  bodies  (A^ . ~ 0 i ^ j)  we  have 


C -*■  — C = — — 
13  2 1 2A 
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C C 

C'12  " Cj+C7  = (A11+A22) 


(2.22) 

(2.23) 


where  Cj  and  C7  are  the  capacitances  of  bodies  1 and  2 with  respect  to 
infinity.  Substituting  Equations  2.22  and  2.23  into  2.20  and  2.21  we  have 
the  usual  result,  namely 


(2.24) 

(2.25) 


Since  A[7  and  A^  in  Equations  2.14  and  2.15  are  positive  (and  because  of 
the  proximity  we  would  expect  Ap  > A} ,)  to  first  order,  the  frequencies 
calculated  from  (2.24)  and  (2.25)  should  be  too  high.  This  latter  statement 
is  true  if  the  electrostatic  treatment  is  the  major  source  of  error  in 
calculating  the  resonant  frequencies.  It  should  be  noted  that  the  frequencies 
described  by  Equations  2.20  and  2.21  did  not  depend  upon  the  particular  shape 
of  the  bodies  but  only  upon  their  mutual  capacitances.  Reference  2 provides 
an  instructive  analysis  for  a particular  system  having  a shape  similar  to 
that  depicted  in  Figure  3. 
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2.3 


CRUDE  ESTIMATES  OF  RESMOD  FREQUENCIES 


In  this  section  we  estimate  the  symmetric  and  asymmetric  frequencies 
expected  from  the  RESMOD  configuration  depicted  in  Figure  1.  The  dependence  of 
the  frequencies  on  general  parameters  such  as  rod  length  and  diameter  should  be 
clear  from  the  way  in  which  the  estimate  is  made.  We  first  estimate  the  capaci- 
tances by  equating  the  area  of  the  two  bodies  to  an  equivalent  sphere.  The  area 
of  the  end  bodies  is  (2) (10.5) (13)+(2) (13) (4.5)+(2) (10.5) (4.5)  or  485  square 
inches.  The  capacitance  Cj  is  then  the  radius  of  the  sphere  or 


area  of  body  . 


in  square  inches  * conversion  to  square  cent imetersj  *^2 


[(485)  T 1/2 

4tJ  (6.45)J  = 15.8  cm  = 17.5  x lo'12  farad  . (2.26) 

The  area  of  the  center  body  is  6(13)2  or  1014  square  inches  so  that 

C2  = 22.8  cm  = 25.3  x io'12  farad  . (2.27) 

The  inductance  of  the  rod  is  given  by 

L = 2 x length  of  rod  in  cm  x £n  (ratio  of  length  to  radius  of  rod)  n.h. 


L = (2) (2.54) (9.5)  x lo'9Hn  C2-!4) (9.5)  = _ 29g  x 1Q-6  h _ (2>2g) 


Given  these  values  of  capacitance  and  inductance  the  corresponding  period  of 
the  symmetric  mode  T . is,  from  Equation  2.24, 

I i r r \ 1/2 


/ LC1C2  V 

VC2*2C1  / 


= 9.3 


(2.29) 


= 6.8  x 10  rad/sec 
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(2.30) 


The  period  of  the  asymmetric  mode  in  is,  from  Equation  2.25, 

,1/2 


and 


t 2 = 2tt(LC1)  ' = 14  n.s. 


w2  = 4.4  * 10  rad/sec 


2.4 


ENERGY  OF  RESMOD  SYSTEM 


In  order  to  obtain  damping  constants  it  is  necessary  to  obtain  the 

i th 

time  averaged  energy  of  the  system.  The  energy  e for  the  £ mode  of  an  N 
body  system  is  given  by 

L _N_  on  1 »— . 0 ~t 

(2.31) 


e*  = i 


1 v*  V"  1 V*,  fT£  -.2 


.£  • 


where  Q is  related  to  the  currents  by  Equation  2.4.  For  the  RESMOD  system, 
where  the  approximations  expressed  by  (2.22)  and  (2.23)  are  valid  we  have 

l £ 2 £21  (Q2)2  1 2 2 

e = [(Qj)  + (Q3)  ] 2T~  + ~2C^~  + \ L^1  + I2)  • (2-32) 

£ 

The  currents  have  the  form  I sintco^+b^) , where  b^  are  phase  constants. 

With  this  form  for  the  currents  it  is  easy  to  show  that  the  energy  averaged 

£ 

over  a cycle  <£  > is 


£ £ 2 

<e  > = L(I*)Z  , i = 1,2  , 


(2.33) 


£ t. 

where  IQ  is  the  magnitude  of  the  current  flowing  in  system  for  the  £in  mode. 
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SECTION  3 

RADIATION  DAMPING 


If  the  energy  radiated  by  a system  is  small  enough  so  that  the 
amplitude  of  oscillation  is  not  changed  appreciably  in  one  cycle,  we  can 
approximate  the  radiation  loss  and  then  approximate  the  change  in  amplitude. 
The  process  of  approximation  is  as  follows:  we  approximate  the  energy  loss 
by  assuming  the  geometrical  distribution  of  current  on  the  radiating  system 
is  known  and  has  the  same  form  as  the  non-radiating  system  (the  same  system 
treated  as  though  it  were  not  radiating) . The  energy  of  the  non-radiating 
system  is  calculated  in  terms  of  its  current  amplitudes  and  then  the  time 
derivative  of  the  energy  of  the  system,  averaged  over  one  cycle,  is  equated 
to  the  energy  loss  also  averaged  over  one  cycle.  The  time  dependence  of  the 
current  amplitudes  is  determined  by  solving  the  resulting  differential  equa- 
tion. We  require  only  that  the  wavelength  of  the  radiation  A be  large  enough 
so  that  our  judgement  about  the  form  of  the  current  densities  is  correct. 

In  effect  this  usually  means  A is  not  smaller  than  A . where 

mm 


2t\V /\ 


min 


(3.1) 


and  £'  is  an  effective  radius  of  the  system  (2ttS, ’ is  the  perimeter  of  the 
system) . 


The  bodies  or  systems  we  are  discussing  are  those  whose  lowest 
modes  have  a wavelength  A larger  than  the  perimeter  of  the  system.  The 
frequency  of  this  mode  can  be  approximated  by  means  of  the  lumped  parameter 
electrical  models  of  the  system  discussed  earlier.  The  form  of  the  currents 
on  the  system  used  in  approximating  the  radiation  damping  depends  upon  this 
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lumped  parameter  quasi-static  model.  For  example,  currents  on  rods  or 
wires  are  assumed  constant  along  the  rods  or  wires  in  a lumped  parameter 
model.  The  lumped  parameter  assumptions  would  not  be  correct  for  a 
higher  mode  where  currents  in  wires  might  not  be  uniform,  for  example. 

As  stated  in  the  previous  paragraph,  the  calculations  in  this  section  will 
be  accurate  so  long  as  X is  not  so  small  that  our  simple  assumptions  about 
the  current  distributions  are  incorrect. 


We  will  be  concerned  in  this  section  mainly  with  calculating  the 
damping  for  the  model  depicted  in  Figure  3.  The  method  can  be  applied  to 
more  general  bodies  and  is  based  on  three  circumstances.  In  summary  these 
circumstances  are:  (1)  the  knowledge  that  the  radiation  rate  is  small; 

(2)  a knowledge  of  the  system  energy  associated  with  the  particular 
frequency  of  radiation;  and  (3)  the  knowledge  of  the  distribution  of  currents. 
The  calculations  to  follow  are  an  example  of  these  circumstances. 

Circumstance  1 is  checked  after  the  calculation  and  serves  as  a consistency 
check  on  to  the  calculation  for  any  particular  set  of  system  parameters. 


It  is  well  known  (a  short  derivation  is  given  in  Appendix  I for 
convenience)  that  the  average  time  rate  of  radiation  of  energy,  <de  /dt>, 

CO 

for  a system  oscillating  with  a particular  angular  frequency,  to,  is 


7 

C / |i)  \ /*  i . 9 a 7 

8?  (c)  J IpJ  U-(P(0-rDsineded<!>  , 


(3.2) 


where 


■*/ 


1 (r')e 

to  ' 


-ir-r'  oo/c 


dr 1 


(3.3) 


18 


I 


In  Equations  3.2  and  3.3,  r,  0,  and  <J>  refer  to  the  coordinates  of  the  field 

A 

point  in  spherical  coordinates,  r is  the  unit  vector  to  the  field  point,  r' 

-y 

is  the  source  point  position,  is  the  source  current  associated  with  the 
angular  frequency  w and  p^  is  the  unit  vector  in  the  p^  direction.  We  first 
calculate  p for  the  model  shown  in  Figure  3 and  then  find  the  rate  of  energy 
dissipation.  To  calculate  p we  first  need  the  currents  on  the  various 
sections  of  the  system.  We  drop  the  subscript  w on  quantities  for  simplicity 
at  this  point.  On  the  two  rods,  labeled  4 and  5 in  Figure  3,  the  current  is 
assumed  constant  along  the  length  and  cross  section. 

3.1  ASYMMETRIC  MODE 

The  first  part  of  our  computation  will  assume  that  the  currents 
in  both  booms  are  in  the  same  direction.  Then 

J4  = J5  = I/ud“ k , «2  < | z | < g + , (3.4) 

where  I is  the  current  in  the  wire  and  d is  its  radius;  the  subscripts  4 
and  5 on  J in  Equation  3.4  refer  to  the  sections  of  the  body  labeled  by  those 
numbers  in  Figure  3.  For  bodies  1 and  3 we  assume  the  current  is  along  the 

A /A 

surface  (in  the  0",  0'  direction)  and  goes  to  zero  at  0"  = - tt  and  0'  = 0 
respectively.  If  A > 2tt«^  and  >>  d the  charge  distribution  is  quasi- 

static and  the  charge  density  is  constant  over  the  sphere.  The  surface 
current  density  J0n(n=l,3)  on  the  spheres  can  be  solved  for  by  means  of  the 
continuity  equation 


1 9 „ • ns  I 

ft  sin0  30  fJ0nSln0)  = 2 

1 4tt  n 


, n = 1 , 3 


and  the  condition  that  the  current  (2nJQ  «.sin0)  goes  to  zero  at  the  proper 

un  1 

angles.  These  solutions  are 


J ! cote^n  g„ 

1 4tt  ft,  ’ 


(3.6) 


x tan8'/2  g, 


4 ir « 


(3.7) 


The  surface  current  on  the  center  sphere  is  found  by  the  condition  that  no 
charge  is  built  up  on  the  surface  so  that 


-1 


J-,  = - (2Tr«2sin0)  10  • 


(3.8) 


Substituting  Equations  3.4  and  3.6  through  3.8  into  Equation  3.3  we  find 
that 

5 


P = Z Pn 
n=l 


(3.9) 


where 


- = 2 e+iuycRcos0  £ j cos30I/2  sin0/2de 

1 C In 


£ — Ie 
c 


iu)/cRcos0 


k , 


and 


-*■  1 T -ico/cRcos0T 

p = — Ie  k . 

r 3 c 


(3.10) 


(3.11) 


Consistent  with  our  quasi-static  assumption  about  the  charge  distribution  on 
the  sphere  we  have  assumed  that  (w/c^)”  < 1 in  Equations  3.10  and  3.11.  We 
have  also  used  the  definition 

R = <■/ + g . (3.12) 
Using  the  assumption  u)/c«9  < 1 for  the  central  sphere  we  have 


p,  c 21  — k 
r3  c 


For  the  rods  we  find,  assuming  wd/c  < 1,  that 


(3.13) 
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p4  = I k fe~ icos8zw/cc 


-(g  + «->) 


Likewise 


= il (ojcosQ) _1j_e+icos0 (u/ c)(&+" 2^  +eicos0(a)/c)  C»2)]{(3. 14) 


P,  = ilCucose)-1^-10086^^^)^-10050^)^]^  (3.i5) 


Substituting  Equations  5.10,  3.11,  3.13,  5.14  and  3.15  into  Equation  3.9 


-+  4 "i  _i 

P - Ik.2  — cos (u)/cRcos@)  +2 (cosOto)  sin(io/c(g+«^)cos0) 


+ 2<t^l c - 2(wcos9)  ^sin^/c^cosO)  j , 


(3.16) 


cos  (cj/cRcos9) 


-2  1 

+ 4(wcos0)  sin^(w/c(g+«?)cos0) 

"1-1  I 

+ 8 — (tocos©)  cos(to/cRcos0)sin(w/c(g+«?)cos0)  | , (3.17) 

where  in  finding  from  p we  have  used  the  approximation  that  u.'«,/c  < 1 
in  the  last  term  of  Equation  3.17  (that  is  we  drop  terms  higher  than  first 
order  in  w^w/c).  IQ  is  the  magnitude  of  I . To  calculate  the  rate  of  energy 
loss  we  need  to  compute  the  integral 


• I 

= / |p|2(l-(P’?)2)sin0d0 


y | p | ^sin30d0  , 


(3.18) 
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if  p only  has  a component  in  the  z direction.  Letting 
x = COS0  , 

then 

dx  = - sin0d0  , 

and 


sin20  = 1 - x“ 


Y(w)  can  be  expressed  as  three  integrals 


Y(w)  = Bj  + B2  + 63  > 


where 


2 1 

= 8("^)  f (l-x^)cos2  (c(jX)dx 
0 


B2  = 8(w)'2 J (x'2-l)sin2(a2x)dx  , 


0 


1 


1 - 1 r -1 

3 = 16  — u>  J (x  -x)cos(a1x)sin(a,x)dx  , (3.25) 

0 


and 


t*1  = (w/c)R  , 


a2  - (u)/c)(g+«2)  = (u)/c)(R-«  ) , 


(3.19) 


(3.20) 


(3.21) 


(3.22) 


(3.23) 


(3.24) 


(3.26) 

(3.27) 


where  we  have  used  Equation  3.12  to  obtain  the  second  form  of  Equation  3.27.  It 
can  easily  be  shown  that 


6:  = jy  + j ocj  (-2a1cos2a]+sin2a1)J  , 


(3.28) 
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r 


f?2  = 8o)  2J-  1 + j [cos(2a?)  + y aj1  sin(2a2) 

+ 2a2Si(2a2)]j  , (3.29) 

and 

a \ _ 1 i 

6-  = 8 — to  j Si  (Oj  +a2)  + Sifo^-ctj) 

+ (oij+ap  1cos(a1+a2)  + (o^-a^^costo^-a  ) 

-2  -2  I 

- (a^a,)  sin(a1+a2)  - (c^-ap  sinf^-a^  . (3.30) 

Equations  3.26  through  3.30  together  with  Equations  3.13  and  3.22  are  those 
necessary  for  calculating  the  radiated  energy.  The  radiated  energy  is 
obtained  by  substituting  Equation  3.18  into  Equation  3.2: 


or 


where 


2,2  , , 

<*>  I0Y(w)  , 


9 


(3.31) 


(3.32) 


c 2 

•tf(<*>)  = ~ 2 — a)  Y(w)  (c.g.s.) 

15ci)2Y(cu)  ohms  . (3.33) 

If  we  def ine  ,y?(oj)  as  the  radiation  resistance  of  the  asymmetric  mode,  the 
normalized  radiation  resistance , ;>p(u>)  is  given  by  the  relation 

.#f(u))  = SOfnj+n^]  2,^P(to)  ohms  . (3.34) 

In  Equation  3.34  n^  is  defined  using  the  radius  of  the  end  body,  , 
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n 


1 


u> 


1 c 


(3.35) 


and  n7  is  defined  using  the  radius  of  the  center  body,  a and  the  length  of  the 
strut,  g, 

n2  s (g+"2)  c • (3’36) 

Plots  of;*P(u>)  appear  in  Figure  4.  The  curves  A to  I correspond  to  various 
values  of  nj.  In  the  low-frequency  limit  (nj-m^OJ.rf+l  and  the  radiation  is 
purely  dipole.  Substituting  from  Equation  2.33  into  Equation  3.32,  where  we 
equate  d/dt<e>  with  <dc/dt>-  wc  find  that 

. •a’M  t 

(lj(t))2  = ( Iq (0) ) “ e 2L  , (3.37) 


or  the  decay  constant  0,  for  the  asymmetric  modal  current  (superscript  1)  is 


a 


4L 


(3.38) 


where  L is  the  inductance  of  one  rod. 


3.2  TESTS  OF  THE  ASYMMETRIC  MODE  DECAY  CONSTANT 


If  we  let  <•/ 1 be  equal  to  zero  we  have  a system  in  which  a rod  is 
connecting  two  spheres.  References  1 and  5 show  how  a rod  loaded  at  both 
ends  by  capacitive  bodies  may  be  treated  as  a rod  with  an  equivalent  length. 
The  decay  constant  formula  for  the  equivalent  length  formulation  is  tested 
in  Reference  5 using  two  computer  simulations  (therein  called  example  2 
and  example  3).  These  simulations  model  the  geometry  depicted  in  Figure  5, 
two  cylindrical  cans  connected  by  a strut.  Figure  5 shows  the  parameters  of 
the  two  simulations  together  with  the  decay  constant  and  frequency 
for  each  simulation.  With  w taken  from  the  computer  simulations  we  will 
test  Equation  3.38  against  og.  To  do  so  we  must  first  obtain  an  equivalent 
spherical  radius  for  the  cylindrical  cans.  Setting  the  area  of  a sphere  equal 
to  that  of  the  cylinder  in  Figure  5 we  find  the  radius  of  the  sphere,  , 
equal  to 
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A-»-n-j  = 10  ^ E -►  n i = .4 

B n-j  = .1  F -►  n i = .5 

C -*■  n-|  = .2  G ^ n]  = .6 

D + n1  = .3  H -c  n-j  = .7 

I -*•  n1  = .8 

Figure  4.  Normalized  radiation  resistance  for  asymmetric  mode 
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(3.39) 


Using  Equation  3.39  together  with  Equation  3.35,  Equation  3.36  (with  </.  = 0) , 
Equation  3.38  and  Figure  4 we  find  that: 

Example  2 

,y7(io)  = 62  ohms  o = 2.9  x 10^  sec 
Example  3 

,*j>(u>)  = 15  ohms  o = 7.2  x 10b  sec  1 

the  damping  rates  are  15  percent  smaller  than  the  code  value  for  example  2 
and  6 percent  larger  than  the  code  value  for  example  3.  Considering  the  nature 
of  the  approximations  which  led  to  Equation  5.39  the  agreement  is  quite  good. 

3.3  DECAY  RATE  FOR  THE  RESMOD  ASYMMETRIC  MODE 

In  Section  2 we  made  a theoretical  estimate  of  the  asymmetric  mode 
RESMOD  frequency.  Ke  now  estimate  the  damping  for  that  mode.  Using  the 
same  procedure  as  used  in  the  examples  of  the  previous  section  we  find 
that,  iij  = .23,  n,  *=  .69,  .#F(a))  = 64  ohm,  a = 5.4  * 107  sec  * and  the  time 
required  to  damp  one  cycle  is  19  n.s.  In  other  words  the  Q for  the  asym- 
metric mode  is  1.4.  The  radiation  rate  is  large  and  the  radiation  resistance 
could  be  somewhat  in  error  due  to  violation  of  condition  (1),  discussed  in  the 
introduction  to  Section  3. 


SYMMETRIC  MODE  DAMPING 


The  method  for  calculating  the  damping  for  the  symmetric  mode  will 
be  very  similar  to  the  calculation  for  the  asymmetric  mode.  The  difference 
is  that  the  currents  in  the  two  struts  are  oppositely  directed  and  the  current 
is  equal  to  zero  along  the  midplane  of  the  central  sphere.  Proceeding  in 
a manner  analogous  to  the  calculations  of  Section  3.2  we  find  that 


" J A 

Pj  + P3  = i2  — I sin((nj+n-,)cos0)k  , 


(3.40) 
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-+  -v  i 4 1 2 n2  ~ 

P4  + P5  = cos9w  Sin  (c°s6T)k  * 


(3.41) 


and 


p2  = 0 , (3.42) 

where  I is  the  current  flowing  in  the  system  and  we  have  used  the  defini- 
tions expressed  by  Equations  3.35  and  3.36  for  n.  and  n.,.  To  find  the 

y 1 7 ^ 

energy  loss  we  must  compute  the  integral  of  |p|  sin  0 from  9 = 0 to  0 = i 

, , 2 i 5+2  2 

(where  |p|“  = |(  l p,  j ) . Defining  this  integral  by  I Y'(od)  (In  is  the 
i = l u u 

magnitude  of  the  current  flowing  in  the  system)  we  find  that 

W2Y’(W)  = n2(frf2)  + 4(f3-f4)  + 4nx(f5-f6)  , (3.43) 


where 


1 

r 2 1 1 

fj  = J sin“(otx)dx  = (a-  — sin(2a))  , a = nj  + n^ 

0 

?2  -J~  x“sin^ (otx)dx  e t ~ (2acos(2a)  + (2oT-l)sin(2a)) 

0 8a"5 

, _ r'sin4  3 n /cos(2„  ) Si(2n  ) ^ 

f3  - J r ’ • 8 - TV  -.i;  2 > 


(3.44) 


n?  ,cos(n  ) , 

T (-T—  + siM  • 


/xn. 


■\dx  = J_  /A  n 

sin(n2) 

sin(2n,)  v 

rX  n2  \ 16  n2 

4 

32  / 

2 (n2x) 

1 „ 

1 , 

sin  — ^ — dx  = 

jS  (a)  - 

j (Si(n2+a)  + 

and 
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= ^'xsin(ax)  sin“|— ^jdx  = ~ -y  (sina-  cosa) 
0 a 


- i (n2+a) 

“ (sin  (n_,+a) 

- (n?+ot) cos (n ,+a) ) 

i , 

-2,  . , 

- (a-n?)cos(a-n,)) 

- j (a-n2) 

(sin(a-n0) 

The  rate  of  change  of  the  average  energy  in  the  system  is  related  to  u>  Y’ (w) 
through  an  equation  analogous  to  Equation  3.32.  In  the  case  of  the  symmetric 
mode  the  radiation  resistance  ,y?'(u))  is 

.tf'  (w)  >=  15uTY' (w)  ohms  . (3.45) 

We  again  define  a normalized  resistance,^' (w)  by  the  equation 

.'/t'  (“)  c SOfnj+n-,)^'  (to)  , (3.46) 

where.^f' (to)  is  plotted  in  Figure  6 as  a function  of  Oj  and  n?.  In  the  low- 
frequency  1 imit  ,J/t'  (to)  -*■  0.  This  is  so  because  the  radiation  from  the 
symmetric  mode  essentially  arises  from  two  oppositely  directed  dipoles;  in 
the  long  wavelength  limit  these  dipoles  are  superimposed. 


2s  2 


Because  the  averaged  energy  in  the  symmetric  mode  is  L(Ip) 


from 


Equation  2.33,  the  decay  constant,  o',  for  the  current  is 
a'  = ' 


. (oo) 

4L 


(3.47) 


except  for  the  primes  this  formula  is  exactly  the  same  as  that  of  the  asym- 
metric mode. 


3.5  DECAY  RATE  FOR  THE  RESMOD  SYMMETRIC  MODE 


As  an  example  of  symmetric  mode  damping  we  calculate  the  decay 

rate  of  the  RESMOD  system.  Using  the  theoretical  estimate  of  the  frequency 

given  by  Equation  2.29  and  the  equivalent  spherical  radii  given  by 

Equations  2.26  and  2.27  we  find  that  n = .36,  n,  = 1 . 1 , ,#>'  (to)  = 16  ohms, 
7-1  1 1 

o'  = 1.3  x 10  sec  and  the  time  required  to  damp  one  cycle  is  77  n.s.  The 
Q for  the  symmetric  mode  is  then  8.2  . 
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1 .EC 


A -*•  n-j  = 10‘3 

B •*  n1  = .1 

C •+  n1  = .2 

0 -*•  n1  = .3 


E n1  = .4 
F - n]  = .5 

G -<■  n1  = .6 

H -*•  n-j  = .7 

I -*■  n1  = .8 


Figure  6.  Normalized  radiation  resistance  for 
symmetric  mode. 


SECTION  4 
SUMMARY 


For  general  systems  whose  lowest  resonant  frequencies  can  be  cal- 
culated by  means  of  electrostatic  and  magnetostatic  concepts  we  have  shown, 
in  the  second  section  of  this  report,  how  the  equations  for  these  frequencies 
can  be  arrived  at  without  the  need  for  an  equivalent  circuit.  (For  those 
not  versed  in  drawing  equivalent  circuits,  a system  like  that  depicted  in 
Figure  2 might  appear  formidable.)  The  method  used  is  useful  for  estimating 
the  errors  made  in  approximating  these  resonant  frequencies.  The  magnetostatic 
parts  of  the  system  were  assumed  not  to  interact  in  the  derivation  of  the 
equations  of  Section  2.  It  would  be  a simple  matter  to  include  these  inter- 
actions in  the  equations,  when  they  are  important. 

Using  a generally  applicable  method  the  damping  rates  of  a symmetric 
three-capacitive  two-inductive  body  system  are  calculated  for  the  lowest  two 
modes,  in  Section  3.  It  is  assumed  that  the  modes  are  actually  normal  modes 
of  the  system  so  that  when  both  modes  are  stimulated  the  radiation  from  one 
does  not  effect  the  amplitude  of  the  other.  The  three  capacitive  bodies  are 
spheres  and  the  whole  system  loosely  corresponds  to  the  RESMOD  geometry  used 
in  the  MRC  Phase  IV  A and  B exploding  wire  experiments.  The  radiation 
resistances  corresponding  to  the  two  modes  are  normalized  to  the  low-frequency 
asymmetric  mode  limit  and  plots  of  the  normalized  resistance  appear  in 
Figures  4 and  6.  These  figures  together  with  Equations  3.34  and  3.46  and  the 
definition  expressed  by  Equations  3.35  and  3.36  can  be  used  to  calculate  the 
actual  radiation  resistances. 
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Numerical  examples  are  provided  in  Sections  2 and  3.  These  examples 

use,  for  the  most  part,  the  actual  dimensions  of  RESMOD  with  rods  of  radius 

.05  cm.  The  numerical  estimates  for  frequency  made  in  these  examples  are 

crude  in  the  sense  that  the  actual  elastances  A. . were  not  known.  A zeroth 

ij 

order  approximation  is  made  where  the  interaction  between  bodies  is  assumed 
zero.  Such  an  approximation  is  valid  when  the  centers  of  the  capacitive 
bodies  are  separated  by  about  3 body  radii.  This  condition  is  just  about 
satisfied  for  RESMOD  but  would  not  be  for  general  systems.  To  apply  the 
methods  outlined  in  Section  2 to  general  systems  an  experimental  and/or 
analytic  procedure  needs  to  be  developed  to  measure  or  compute  the  elastances. 
These  procedures  could  take  the  form  of  electrical  tests  or  moment  method 
computer  codes. 
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APPENDIX  I 

RADIATION  FROM  A SYSTEM  OSCILLATING  WITH  ANGULAR  FREQUENCY 


In  this  appendix  we  integrate  the  Porting  vector  over  a sphere 
at  infinity  and  then  average  over  one  cycle  to  obtain  the  average  rate  of 
energy  loss  due  to  radiation.  The  Poynting  vector  is  defined  as 


where 


and 


? c + -+ 

4?  E * B ’ 


I M 

c 8t 


B = V x A . 
In  the  Lorentz  gauge 


Cl.  l) 


(1-2) 


(1-3) 


and 


1 

CJ  \t-P\ 


dr' 


(1.4) 


If 


<P  Cr 


,t)  = / <l_ 

J i?  - ?• i 


dr ' 


^(?,t)  = RP?n(?')eiwt  , 


(1.5) 


(1.6) 
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1 


where  Rp  denotes  the  real  part  of,  then  for  large  r 

iu(f)  J 


A = Rp 


cr 


where 


- f t -iw/c  r-r'  ■* 
p = I JQe  dr'  , 


and 


We  also  have 


t ' = t - r/c  . 


iwt ' 


= Rp 


where 


P 


„ _ , -ir-r'  co/c 
q = / pQe  dr'  . 


noting  that 


iwq 


- - flV'.  “/c  dr' 


(1.7) 


(1.8) 


(1-9) 


(I. 10) 


(I. ID 


(1.12) 


from  the  continuity  equation,  we  have  after  integrating  Equation  1.12  by 
parts 


q = 


(1.13) 


Substituting  Equation  1.13  into  I. 10  and  the  result  together  with  Equation 
1.7  into  Equation  1.2  we  have 

icot ' iojt ' 

E = Rpi  ^ (p_r.£)  = Rpi  | S__  r x (jjx?) 


w 1 , . iwt 


7 (ie"*'*’  r x (pxf)  . ie”1U)t ' r x (p*xr)).  (1.14) 


36 


From  substituting  Equation  1.7  into  Equation  1.3  we  also  obtain 


d o • w e 
B = Rpi  - — 


2c  F (ie  PXr  - ie 


-iwt ' 


p**r) . 


(1.15) 


Substituting  Equations  1.14  and  1.15  into  Equation  1.1,  integrating  over  the 
sphere  at  infinity  and  averaging  over  one  cycle  we  obtain  the  average  rate  of 
radiated  energy  <d£  /dt>  at  the  frequency  u>: 


8tt 


(pxr) ) x 


(Pxr))  • r dft 


* ~ Wirf  IpI~(?xP)  • (*XP)  dft 

2 

= " gF  f !p|2(1  ' (P'f)“)  dil  . (1.16) 
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